Effects of randomness on the critical temperature in quasi-two-dimensional organic 

superconductors 
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The effects of non-magnetic disorder on the critical temperature T c of organic weak-linked layered 
superconductors with singlet in-plane pairing are considered. A randomness in the interlayer Joseph- 
son coupling is shown to destroy phase coherence between the layers and T c suppresses smoothly 
in a large extent of the disorder strength. Nevertheless the disorder of arbitrarily high strength 
can not destroy completely the superconducting phase. The obtained quasi-linear decrease of the 
critical temperature with increasing disorder strength is in good agreement with experimental mea- 
surements. 
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I. INTRODUCTION 

Organic molecular crystals n-(BEDT-TTF) 2 X [ab- 
breviated as k — (ET)2X] are in the center of attention 
due to their unusual normal metallic and superconduct- 
ing properties 1 . The flat ET molecules in n — (ET) 2 X 
organic metals dimerize to form molecular units that 
stack in planes on a triangular lattice 2 . The anions 
X, which modify from X = Cu[N{CN) 2 ]Cl through 
X = Cu[N{CN) 2 ]Br to X = Cu{NCS) 2 , separate the 
planes and accept one electron from each BEDT — TTF 
dimer. Most of the ET-based superconductors (SCs) are 
strongly anisotropic quasi-two-dimensional (quasi-2D) 
conductors because the conductivity is approximately 
isotropic within the layers of the ET donor molecules 
but smaller by a factor of ~ 10 3 in the perpendicular 
direction. Measurements of the superconducting coher- 
ence lengths 3 within- £y and perpendicular £j_ to the su- 
perconducting planes in e.g. k — (ET) 2 Cu[N(CN) 2 ]Br 
yield £|| « 37 A and £j_ ss 4 A, the latter of which is much 
smaller than the interlayer distance ~ 15A. This fact 
suggests that superconductivity in the direction perpen- 
dicular to the plane may involve Josephson tunneling. 

Low temperature properties of organic SCs are known 
to be very sensitive to disorder 4 . Alloying with an- 
ions, x-ray irradiation, or cooling rate controlled anion 
reorientation introduces non-magnetic randomness into 
the system, however leaving unchanged, to a large ex- 
tent, in-plane molecular structures. Recently, the ef- 
fects of non-magnetic disorder on superconductivity in 
organic n — (ET) 2 Cu(SCN) 2 have been studied exper- 
imentally in Refs. 5,6 . The non-magnetic disorder was 
introduced in these experiments via irradiation by ei- 
ther x-rays or protons 5,6 and via partial substitution of 
BEDT -TTF molecules with deuterated BEDT -TTF 
or BMDT—TTF molecules 6 . All disorder seems to affect 
the terminal ethylene group and anion bound structures. 
The measurements for samples with molecular substitu- 
tions show 6 that the mean free path I is longer than the 
in-plane coherent length £|| , indicating that the supercon- 
ducting planes can be considered to be in clean limit. T c 



was found 5 to fall quasi-linearly with defect density, and 
the dependence exhibits a sharp change in slope from 
0.31 to 0.15 around a threshold value of the interlayer 
residual resistivity p$ « 2 flcm. The main feature of 
the experiments is that the samples exhibit a supercon- 
ducting ground state even at the highest defect densities, 
and there is not a SC-normal metal phase transition dif- 
ferent from quasi-lD organic SCs , where the random- 
ness transforms the system into a normal metallic state. 
In the light of the experimental data, the Abrikosov- 
Gor'kov's theory 8 for non-magnetic defects in non-s-wave 
SCs seems to fail to explain the experimental data. 

We study in this article the effects of randomness in the 
Josephson coupling energy on the critical temperature 
of weak-linked quasi-2D SCs. Therefore, the influence 
of a possible in-plane molecular disorder on the super- 
conducting properties of the system is ignored. Suppres- 
sion of superconductivity in the presence of non-magnetic 
impurities can in general be realized by destroying ei- 
ther the modulus or the phase coherence of the order 
parameter 9 ' 10 . Although strong fluctuations of the order 
parameter phase destroy off-diagonal long-range order 
(ODLRO) in an isolated superconducting plane 11 , the 
point-like topological defects of a "phase field" such as 
"vortex" and "antivortex" of Kosterlitz and Thouless 12 
sets up a quasi-long range order in the system. 



II. CLASSICAL PHASE FLUCTUATION 
REGIME 

The strongly anisotropic organic SCs with in-plane sin- 
glet pairings are modeled as a regularly placed supercon- 
ducting layers with Josephson-coupling between nearest- 
neighboring layers with a classical free energy functional 
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where <Pj(r) denotes the phase of the order parameter 

(2) 

Aj(r) = \Aj\exp(i(pj(r)), Ns (T) is the surface den- 



sity of superconducting electrons; Ng 2 \T) = N [ ^'(0) 
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The last term in Eq.(l) 



describes the Josephson coupling with the energy i^jj +g 
Fluctuations of the order parameter modulus can be 

(2) 

neglected for pure SCs far from T^>, the mean-field crit- 
ical temperature calculated for an isolated layer. There- 
fore, the contributions to F st {ip} in Eq. (1), coming from 
the modulus of the order parameter |Aj|, are omitted. 



We assume the Josephson energy E 



3,3+9 



to be 



dom parameter with Gaussian distribution, centered at 
the mean value E g , given by 

P{E jtj+g } = (2nW 2 )- 1 / 2 e^{-(E j , j+g -E g f/(2W 2 )}. 

(2) 

W 2 is taken as a measure of disorder strength. 

Employing the replica trick one can calculate the av- 
erage value of the order parameter cos(^j(z)) 
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where Z = J Dip exp (— ^F st {ip, rf\) is the partition func- 
tion with respect to the free energy functional F st {(p, r/} 
which contains, in addition to Eq. (1), the generat- 
ing field term J2j Vj cos {fj( r ))- The double bracket 
((. . . ))dis is a shorthand notation for the double average 
over thermodynamic fluctuations and over disorder. 
Integration out the Gaussian random variables yields 



((cos(^(r)))) dis 
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WN {2) (cos[y jo - <p jo+go ]) 2 - (cos 2 [y 30 - ip jQ+go ]) 
T (cos[^ - ip jo+go }) 
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The expression for the effective free-energy functional 
at the saddle-point Q , go , given by Eq. (6), becomes 
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similar to that for a regular quasi-2D SC with renor- 
malized inter-layer Josephson energy, E g — > E g — 

W 2 Nf> (cos[y 30 -y J0 + g0 ]) 2 -(cos 2 [y_,- -y J0 + g0 ]) 

Equation for T c in quasi-2D SCs is derived from Eq. (4) 
by using the self-consistent mean- field method 13 , which 
consists in replacing the cosine term of Eq. (5) as 

^£ g cos(pj - ipj+g) -> ry£ ± ((cos(»)) c cos </?(», (7) 



where rj is the coordination number, and E± rts t^/ep 
with t± and ep being the interlayer tunneling integral 
and the Fermi energy. The phase correlations on the 
nearest-neighboring layers in this approximation are sim- 
plified by describing them as a motion of a phason in 
the average field of phases with the most probable value, 
which coincides with the average value for a clean sys- 
tem ((cos((p))) c = ((cos(ip))) di S 13 ■ For a dirty system 
the most probable value differs strongly from the aver- 
age value. Indeed, we assume that a distribution func- 
tion of the order parameter in the presence of random- 
ness is broad and asymmetric. This broadness and asym- 
metry becomes stronger around the critical temperature 
due to huge thermal fluctuations. Therefore, knowledge 
of the arithmetic average is insufficient, and infinitely 
many moments give a contribution to the distribution 
function of the order parameter at the tail. We identify 
{(cos((f))) c with the most probable or typical value of 
the order parameter. For the disordered SC we choose 

«cos„» c = «cos^ s - <<c ° sy) ^t Sy> ^" which 
resembles a change made by the saddle-point (6) in the 
free-energy functional. The functional integral over the 
phases in Eq. (4) can not be evaluated yet, even after 
this simplification. Taking advantage of the smallncss of 
E±((cos(ip))) c near T c , however, an expansion of the in- 
tegrand of Eq. (4) in this quantity allows us to obtain the 
following equations for ((cos(</?))) dis and ((cos(<^)) 2 ) dis 

riNs E± f 

({coa(tp)))dis = fc S T \ d 2 r(cos ip(0) cos ip(r)) (cos ip) c 
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((cos(^)) 2 ) dis = (1 + W 2 /El ) ((cos ip)) dis . (9) 

The final equation for T c , obtained from the above writ- 
ten expressions, reads 



where Q g (r) denotes a Hubbard-Stratonovich auxiliary 1 
decoupling field. In order to clarify a character of the 
saddle-point for the variable Cj,g( r )i one writes the ex- 
pression (4) for ((cos[(fij(r) — ipj +g (r)])) d is and find the 
saddle-point Cj ,ga as 
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The phase-phase correlator in Eq. (10) is calculated in 
the clean limit of the 2D free energy functional, obtained 
from Eq. (1) by setting £j j+ g = 0, which yields, 9,11 
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where £|| = 2 h7VF (2) with In 7 = c = 0.577 is the in-plane 

coherence length. Real-space integration of the corre- 
lator (11) in Eq. (10) for the critical temperature im- 
poses the following restriction on the critical temperature 
n ik f; T i2 h + 2 < yielding T C >T*, where 



l/T* = + 2k B /e F . 



(12) 



T* may be identified as the Kosterlitz-Thouless transi- 
tion temperature. The order parameter's phases correla- 
tion between nearest-neighbor layers disappears as T c ap- 
proaches T*, and system reveals effectively 2D supercon- 
ducting behavior at T = T* . So, the critical temperature 
in quasi-2D SC varies in the interval of T* < T c < T c (2) . 

One can estimate T* for k - (ET) 2 Cu(SCN) 2 SC with 
(2) 

Tc ~ 10. 5K for which the Fermi velocity and the ef- 
fective mass of electrons are measured 5 ' 6 to be Vp ps 
4 x 10 4 m/s and m* ps 3 mo, respectively, where mo is a 
free electron mass. These data yield T* ss 8.8K, which 
agrees well with maximally dropped critical temperature 
with disorder in Refs. 5,6 . The randomness in the Joseph- 
son energy in the presence of the order parameter phase 
fluctuations destroys the transverse stiffness in the sys- 
tem. Equations (10) and (11) yield 



1 = 



qt 2 {4(1 - e-V«) + 1/(1 - t )}, (13) 



where t is a dimensionless T c -shift, < t < 1, introduced 

as 



t = (e F /2k B ) (l/T c - 1/T C < 
and q is a dimensionless parameter 
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with In 7 = c = 0.577. q decreases from its maximal value 

q = qo = -^S- ( k ^(2) ) ^° zer0 as ^ ne disorder param- 
eter x — W 2 /Ej_ increases from zero up to the maximal 
value x = 1 for strong randomness W ~ The nu- 
meric solution of Eqs. (13) and (14) for the dependence 
of T c on x is depicted in Fig. 1. 

Equation (13) is solved in two asymptotic limits. For 
< t < 1/4, which corresponds to a weak disorder limit 

(2) 

when T c varies around T c , the exponential term is ne- 
glected, yielding 



1 
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In the limit 1/4 < t < 1, which corresponds to rela- 
tively strong disorder limit when T c varies around T*, 
the exponential term in Eq. (13) is expanded yielding 
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FIG. 1: (Color online) Dependence of T c on a: = W 2 /E± 
for go = 0.6, 1.0 and 2.0 is depicted by dashed (green), dot- 
dashed (blue) and solid (red) curves, correspondingly. T c (q) 
increases with qo (or E±), and reaches T* at highest value of 
the randomness W = E± when x = 1. Insert shows T c (x) 
dependence according to the Abrikosov-Gor'kov's digamma 
function, which vanishes for higher values of x, instead of 
saturation in our case. 



III. QUANTUM PHASE FLUCTUATIONS 
REGIME 



The results, obtained above for the classical fluctua- 
tions are valid for weak randomness, when T c and E± 
have relatively high values. Quantum fluctuations have 
to be taken into account for a strong disorder (small T c 
and E±) limit. The effects of randomness in the pres- 
ence of quantum phase fluctuations can be studied by 
starting from a Hamiltonian of weak-linked metallic lay- 
ers with in-layer attractive electron-electron interactions. 
Integrating out the electronic degrees of freedom in the 
partition function, by following the method of Ambe- 
gaokar et al. 14 , yields the following the expression for 
the dynamical free energy functional 



h 

8V 



i,3 J J 
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where the phases depend now on the imaginary "time" 
t too, the dot on the phase means a "time" -derivative, 
Kij is the susceptibility 9,14 . F^{ip} is the stationary 
part of the dynamical free energy functional, which differs 
from the classical functional (1) by additional integration 
over the imaginary "time" r. Repeating the procedure 
of derivation of Eq. (10) for the critical temperature in 
the case of the classical fluctuations, one arrives at 



l = ^xJV< 2 Wl- J-J x 

x J dr j d 2 r(cos[</?(0,0)]cos[(p(r,T)])o.(19) 
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FIG. 2: (Color online) Dependence of T c on x = W 2 /E 2 ± , 
according to Eq. (22), for a — 0.5 and three different values 
of go: go = 0.6 dashed (green) curve, go = 1.0 dot-dashed 
(blue) curve and qo = 2.0 solid (red) curve. The dependence 
of T c on the residual resistivity po in the insert is plotted 
by combining T c (po) suppression in the classical fluctuations 
regime at po < Po — 2Qcm with a slow suppression of T c (po) 
in the quantum phase fluctuations regime (corresponding to 
dot-dashed curves with go = 1) at po > Po- 



The phase-phase quantum correlator for a pure 2D SC is 
calculated to yield 

(cos[^(r,T)-p(0,0)]) = 

= eX p y y 2[1 ; 2 TOS(k ' r = "" r)] | (20) 

A straightforward calculation results in 
(cos[^(r,r)-^(0,0)]) = 



exp 
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the dimensionless quantum parameters. Although a 



4trk B T^ 



epT i/2 ~ a o is proportional to the dynamical parameter 



"0 = nZ 



1/2 



(2h' i /m-')K > wn i cn characterizes a quan- 
tum charging effect in the system, the other parameter 



depends inversely on cto , nevertheless the 



product a/3 does not depend on a . The quantum corre- 
lator (21) becomes the classic one (11) for cxq — > (a <C 1, 
P » 1, but a/3 — > const). Equations (19) and (21) yield 
a dependence of T c on E± and W for two limiting cases 
P > 1 and /3 < 1. In both cases the integration of the 
correlator (21) over coordinates imposes the restriction 



T* < T c < T c (2) (or a/3 > 2), where T* is defined by 
expression (12). For /3 < 1, which corresponds to strong 
charging regime, Eq. (19) after integration over r and r 
results in non-linear equation for t 



t 6 



i - * 



exp(-2 v / a/t ) } = 1, (22) 



where, a = — 51 characterizes the dynamic effects 

too, since a « a^. The numeric solution of Eq. (22) 
is given in Fig. 2 for a = 1.2 and different values of 
go ■ The quantum fluctuations reduce T c considerably and 
alter the results of the classical fluctuations regime at low 
temperatures and small E±_ , which corresponds to strong 
randomness or high resistivity in the T c {x) dependence. 
The slope of, e.g. the dashed (green) curve with qo = 0.8 
changes from 0.54 in the interval < x < 0.5 of Fig. 

1 for the classic fluctuations regime to the value of 0.30 
in the interval 0.5 < x < 1 of Fig. 2 for the quantum 
fluctuations regime, the ratio of which (1.8) is comparable 
with that (~ 2) estimated for the experimental curve 5 . 

For j3 > 1 Eqs. (19) and (21) are solved in the limit of 

2 < a/3 < 8 yielding 

T c = T*(l + q)/{l + T*q/T^}. (23) 
The case of f3 > 1 and a/3 > 8 results in 

1 = ^ {4+1/(1-01, (24) 

an approximate solution of which is given by Eq. (16). 
The case of /3 > 1 or a < 1 corresponds to the weak 
quantum fluctuations limit, and, therefore, the results 
do not depend on the dynamical parameter ao- 

A detailed comparison of the results with the ex- 
periments needs to express T c on the interlayer resid- 
ual resistivity po — irh 4: /(2e 2 m*a±t 2 L Tt), 4: where a± is 
the interlayer distance. The inelastic scattering time 
h/r t — irci mp N(o)\U\ 2 is assumed to relate with a mea- 
sure of the randomness x ~ W 2 as W 2 — ira mp \U\ 2 /2. 
It is necessary to take into account that the charging ef- 
fect in the quantum fluctuations regime reduces the value 
of the interlayer tunneling integral t±_ (or the transverse 
rigidity) 9 . Therefore, po in the quantum fluctuations 
regime is rescaled for a given value of x to a larger inter- 
val in comparision with that in the classical fluctuation 
regime. The dependence of T c on p Q is given in the insert 
of Fig. 2. 



IV. CONCLUSIONS 

In this paper we report disorder effects on T c of quasi- 
2D SCs with random Josephson coupling. The interplay 
of non-magnetic disorder with quantum phase fluctua- 
tions becomes a central factor in suppression of the su- 
perconducting phase in organic quasi-2D SCs. A ran- 
domness in the interlayer coupling energy is shown to 



■5 



decrease T c quasi-linearly, nevertheless the superconduct- 
ing phase does not completely vanish even at arbitrarily 
high strength of the disorder. The present theory ex- 
plains very well the recent experimental measurements 
given in Refs. 5 ' 6 . We neglect in this article effects of in- 
plane disorder on T c in organic SCs. Such randomness 
results in suppression of T c due to the Anderson local- 
ization for non-s-wave pairings, and it seems to destroy 
the homogeneity of the order parameter modulus lead- 
ing to the formation of a cluster-like "superconducting 
island" inside the metallic phase. On the other hand 
the in-plane disorder may "pin" the Kosterlitz-Thouless 



topological defects and destroy the quasi-long range order 
in the system. All these effects deserve further investiga- 
tion. 
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